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ABSTRACT 

Context. Once ALMA full polarization capabilities are offered, it will become possible to perform detailed studies of polarized maser 
emission towards star-forming regions and late-type stars, such as (post-) asymptotic giant branch stars and young planetary nebulae. 
To derive the magnetic field orientation from maser linear polarization, a number of conditions involving the rate of stimulated emis- 
sion R, the decay rate of the molecular state T, and the Zeeman frequency gO need to be satisfied. 

Aims. The goal of this work is to investigate if SiO, H 2 and HCN maser emission within the ALMA frequency range can be de- 
tected with observable levels of fractional linear polarization in the regime where the Zeeman frequency is greater than the stimulated 
emission rate. 

Methods. We used a radiative transfer code to calculate the fractional linear polarization as a function of the emerging brightness 
temperature for a number of rotational transition of SiO, H 2 and HCN that have been observed to display maser emission at submil- 
limetre wavelengths. We assumed typical magnetic field strengths measured towards galactic star-forming regions and circumstellar 
envelopes of late-type stars from previous VLBI observations. Since the Lande g-factors have not been reported for the different rota- 
tional transitions we modelled, we performed our calculations assuming conservative values of the Zeeman frequency for the different 
molecular species. 

Results. Setting a lower limit for the Zeeman frequency that still satisfies the criteria gCl > R and gCl > T, we find fractional polariza- 
tion levels of up to 13%, 14% and 19% for the higher J transitions analysed for SiO, H 2 and HCN, respectively, without considering 
anisotropic pumping or any other non-Zeeman effect. These upper limits were calculated assuming a magnetic field oriented perpen- 
dicular to the direction of propagation of the maser radiation. 

Conclusions. According to our results SiO, H 2 and HCN maser emission within the ALMA frequency range can be detected with 
suitable linear polarization to trace the magnetic field structure towards star-forming regions and late-type stars even if the detected 
polarization has been enhanced by non-Zeeman effects. 

Key words, stars: AGB and post-AGB - Masers - stars: Magnetic Field - Polarization - stars: Formation - submillimeter: stars 



1. Introduction 

Polarized maser emission has been detected towards star- 
forming regions (SFR) and expanding circumstellar envelopes 
(CSE) of late-type stars such as (post-) asymptotic gi ant branch 
(AGB stars and young p l anetary nebulae (PNe) (e.g . lAlves et aTl 
l2012tjAmiri et al.ll2012l;lLe7l -Ferreira et al.l l2012l Vlemmings 
et al. I201 lL l2006bh . Both single-dish and interferometric ob- 
servations have revealed that silicon monoxide (SiO), water 
(H2O), hydrogen cyanide (HCN), hydroxyl (OH), and methanol 
(CH3OH), among others, can naturally generate polarized maser 
emission in t hese enviroments (e . g. Vlemming s et al.ll2 006b: 
Surcis et al.l I20TI iFish & Reidl 120071* iHeroin et al.l 12006: 



Kemball et dj2009). In particular, interferometric observations 



of the masers at radio wavelengths have become a useful tool 
for studying t he magnetic field in and around SFR s and late- 
type stars (e.g. lFish & Reidl2007aHAmiri et al.l2010l) . At shorter 
wavelengths, new instruments will enable the study of maser ra- 
diation from higher vibrationally-excited rotational transitions. 
In particular, the Atacama Large Millimeter/submillimeter Array 
(ALMA) has recently started the first scientific observations. 
Soon all of its capabilities, including polarimetry, will be avail- 
able, providing more than an order of magnitude improvement 



in sensitivity and resolution. In the ALMA frequency range, a 
number of SiO, H2O and HCN maser ttansitions that belong to 
vibrationally-excited levels up to v = 3 can be detected. The 
population inversion of rotational states belonging to different 
molecular species, or within different vibrational levels of a par- 
ticular species, depends on the physical conditions of the emit- 
ting regions. By observing the maser polarization at different fre- 
quencies towards an extended source, we can hence constrain 
the magnetic field properties, i.e. field strength and/or direction, 
throughout large regions within the observed source. 

According to the Zeeman interpretation, the maser polar- 
ization strongly depends on the ratio between the Zeeman fre- 
quency (gfi), the rate of stimulated emission (R), and the rate 
of the decay of the molecular state (F). The Zeeman splitting 
induced by a magnetic field depends on the molecule's shell 
structure. For example, since SiO, H2O and HCN are non- 
paramagnetic, closed-shell molecules, gQ. is expected to be less 
than the intrinsic line wid th Aa>. Maser polarization theory (e.g 
IWestern & Watsonll 19841) predicts fractional linear polarization 
(Pl) levels of up to 100% for 7=1-0 rotational ttansitions 
of diatomic molecules. These polarization levels can be reached 
in the presence of a magnetic field of a few Gauss in, for exam- 
ple, the SiO masing region in CSEs of late-type stars. But in the 
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case of higher rotational transitions (i.e. 7 = 2-1, 7 = 3-2, 
etc.) theory predicts that unless anisotropic pumping is involved, 
the fractional linear polarization should decrease as the angular 
momentum number of the involved state increases. 

Before the observed maser polarization can be related to an 
intrinsic magnetic field, it is necessary to evaluate the ratios be- 
tween the maser parameters R, gQ and T for each single rota- 
tional transition detected with angular momentum higher than 
7 — 1. Here, we determine if the polarized maser radiation 
produced by molecular transitions of SiO, H2O and HCN in 
the ALMA frequency range could generate detectable levels of 
fractional linear polarization, while still satisfying the criteria 
(gQ > R and gQ > F) for which the polarization direction is still 
directly related to the magnetic field. These maser transitions can 
then be used to infer the magnetic field structure towards SFRs 
and in the CSEs of (post-) AG B stars. To do this, we have run 
numerical models adapted from Ned oluha & Watsonl (1 19921) to 
calculate the fractional linear polarization level that can be gen- 
erated by the interaction of the molecular states of SiO, H2O and 
HCN in the ALMA frequency range, with a magnetic field in the 
masing region. 

In section 2 we present a brief review of the maser polariza- 
tion theory and in section 3 we describe the different maser tran- 
sitions of SiO, H 2 and HCN that can be detected in the ALMA 
frequency range. The results of our models for each molecular 
species are detailed in section 4. In section 5 we give an estimate 
of the polarized fluxes that can be expected when observing with 
ALMA. 

2. Maser polarization theory 

Polarized maser radiation provides a unique tool for studying 
the role of magnetic fields inside high-density enviroments such 
as star-forming regions and circumstellar envelopes of evolved 
stars. Although the general maser emission mechanims from as- 
tronomical sources h ave been wel l understood throughout the 
past thirty years (e.g. lElitzurlll992l lGravll2012h . the generation 
and radiative transfer of polarized maser emission have been dif- 
ficult to deal with. The polarization properties of maser emission 
strongly depend on the radiative conditions of the region where 
the maser is being generated (saturated or unsaturated) and on 
the nature of the molecular species generat ing the maser emis- 
sion (paramagnetic or non-paramagnetic) (Watson 2008; Dinh- 
v-Trung l2009T ~ 

Maser-polarized radiation can be produced in both unsaturated 
and saturated frames. Maser emission is considered to be satu- 
rated when the rate for the stimulated emission R overcomes the 
decay rate of the molecular state involved, T. In this case, the 
growth of the possible polarization modes is determined by the 
population of the molecular state s that can interact with a par - 
ticular mode of the po larization (Nedoluha & Watson 1990a). 
iGoldreich et al.l dl973h first identified two regimes where polar- 
ized emission can be generated in the saturated frame for the 
molecular rotational transition 7 = 1 - 0: a) The strong magnetic 
field strength regime, where the Zeeman frequency gQ. » R 
and gQ. » T; and b) The intermediat e magnetic field strengt h 
regime, where (j?Q) 2 /r » R » gQ.IWestern & Watsonl (fl 984). 
iDeguchi & Watsonl dl990h andlNedoluha & Watsonl (Il990al) ex- 
tended the treatment of IGoldreich et al.l (1 19731) . solving the radia- 
tive transfer equations for polarized radiation as a function of the 
emerging intensity R/F, i.e. as a function of the saturation level. 
They considered linear, non-paramagnetic molecules, including 
the rotational transitions 7 = 2 - 1 and 7 = 3-2. Nedoluha & 
Watson (1990a) have shown that in the presence of a plausible 



magnetic field and without differences in the population of the 
magnetic substates involved in the maser emission, the upper 
limit of fractional linear polarization pi that maser transitions 
other than 7=1-0 can achieve is 33%. The fractional linear 
polarization decreases for transitions with higher angular mo- 
mentum, and is a function of the angle 9 between the magnetic 
field (B) and the direction of the maser radiation (k). 
Molecular states with angular momentum higher than 7=1-0 
reach their maximum pi value when the magnetic field lines are 
perpendicular to the direction of the propagation of the maser ra- 
diation, and have pi = when 9 has the critical value 9 cr ~ 55°, 
also known as the van Vleck angle, and when 9 = 0°. 
The angle <p between the polarization vector and the plane k ■ B 
also varies as a function of the ratio gQ/R. For a fixed 9, the 
vector of polarization should be either perpendicular or parallel 
to the k ■ B plane if the condition gQ » R is satisfied. On the 
other hand, for higher values of R, but still in the intermediate- 
strength magnetic field regime, the vector of polarization is al- 
ready rotated away from the k ■ B plane, i.e. if R > gQ the 
vector of polarization is neither parallel nor perpendicular to the 
plane k ■ B. In this context, the van Vleck angle is the limit for 9 
where the polarization vector changes from being parallel to be 
perpendicular to the k ■ B plane. Therefore, as long as the con- 
ditions gQ > R and gQ > F are satisfied, the information about 
the morphology of the magnetic field can be extracted from the 
polarization vector. 

Observationally, fractional linear po larization of up to 100% 
for 7=1-0 SiO maser transitions (Ami ri et al 1 12012b and 
values higher than 33% for SiO molec ular transitions involving 
higher angular momentum states (e.g IVlemmings et al.l 120 1 lb 
have been detected. Anisotropic pumping seems to be a possi- 
ble pro cess to explain such high pi levels dNedoluha & Watsonl 
1990a). It can be produced by differences in the angular distri- 
bution of the radiation field involved in the population inversion 
process. In the absence of a magnetic field, anisotropic pumping 
could produce highly linearly polarized emission. In contrast to 
the effect of a magnetic field alone in the maser region, the frac- 
tional linear polarization can increase with the angular momen- 
tum due to an anisotropic population of the involved magnetic 
substates. Nevertheless, even considering that the linear polar- 
ization has been enhanced by anisotropic pumping, if gQ > R 
and gQ > F for the detected line, the magnetic field structure can 
still be traced directly from the polarization vect or, because th e 
magnetic field is the dominant axis of symmetry (Watson 2002). 
The study of magnetic fields from maser linear polarization thus 
depends critically on an analysis of the Zeeman frequency gQ 
and the stimulated emission rate R. 



3. Observations of (sub)millimetre masers 

3.1. SiO 

Although it was first d etected towards the Orion molecular cloud 
(Snyder & Buhl 1974), later surveys ha ve shown that SiO m aser 
emission is quite uncommon in SFRs (Zapata et al. 2009])- To 
date, the only submillimetre SiO masers detected in SFRs are 
the v = 1, 7 = 1 - 7 = 3 - 2 and the v = 2, 7 = 1 - 
rotational transitions dBuhl et al.lll974t iDavis et al.lll974l) . In 
contrast, strong SiO emission arises from the innermost regions 
of the CSE of late-type stars, and a number of maser transi- 
tions have been detected in vibrational levels of up to v = 4, 
with rotational transitio ns as high as 7 = 8 - 7 (M enten et al.l 
2006, Humphreys 2007 and references therein). Interferometric 
obervations have revealed that the SiO maser emitting regions 
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Table 1. Spontaneous emission coefficient and frequencies of the detected SiO maser transitions in the ALMA frequency range 
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Table 2. Einstein coefficient and frequencies of the detected H2O and HCN maser transitions at submillimetre wavelengths. 
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Notes. (<I) No spectroscopic data available 



form ring-like structures centred on the star, between the stellar 
photosphere and the dust-forming region ~ 2-6 R+ (Diamond 
et al. ll994l : lBoboltz & Diamondll2005l) . The most common SiO 
maser transitions detected are the two lowest rotational transi- 
tions of the v — 1 vibrational level, peakin g around 43 GHz and 
86 GHz (e.g. lBarvainis & Predmo"relll985l) . The list of SiO rota- 
tional transitions, including the isotopologues 29 SiO and 30 SiO, 
which can be observed as maser emission in the AL MA fre- 
quency range, are listed in Table[TJ(Muller et al. 2001). 
Polarized SiO maser emissi on has been detected tow ards 
oxygen-rich (post-)AGB stars. iBarvainis & Predmoreldl985l) re- 
ported the detection of strong polarized emission in the rota- 
tional transitions / = 1 - and J — 2 - 1 , of both v = 1 and 
v = 2 vibrational levels. 

Fractional linear polarization levels between 15% and 40% are 
commonly detected, though levels approaching 100 % have also 
been reported (e.g. Barvainis & Predmore 1985; Kemball & 
Diamond 1 19971: 1 Amiri et al.ll2012h . lShinnaga et all d2004 and 
IVlem mings et al.l d201 ll) reported the detection of high fractional 



linear polarization levels (pi > 40%) of the 7 = 5-4 rotational 
transitions of the v = 1 vibrational level, towards the supergiant 
VX Sgr. Such high polarization levels cannot only be explained 
by the presence of a large-scale magnetic field permeating the 
masing region, but need to be enhanced by non-Zeeman effects. 



3.2. H 2 

Water is among the most abundant molecules in the envelopes of 
cores embedde d in star-forming regions and in the CS Es of late- 
type stars (e.g. Waters et al. 1980; Menten et al. 1990a; Maercker 
et al. l2008l) . At submillimetre wavelengths a number of molec- 
ular transitions of H2O have been detected displaying high flux 
density values and narrow line shapes characteristic of maser 

J 1 1 I 1 1 1 1 I 

emission lines (Waters et al. 1980; Menten et al. 1990b; Yates 
& Cohen 1996). The most studied water maser emission is the 
low-frequency 6i6-523 transition at 22.2 GHz. Since this transi- 
tion is not affected by the atmospheric precipitable water vapor 
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(PWV), it has become a reference in the study of water maser 
emission from astrophysical sources. It has been detected to- 
wards the expanding CSEs of late-type stars and high-velocity 
outflows generated in the envelope of both protostellar objects 
and (post-) AGB stars. This variety of scenarios gives us an idea 
of the broad excitation conditions of the 22.2 GHz maser (high- 
density rej*i£nsji£ 2 _>Ji)ii^rr^ K > 7 > 
200 K. lNeufeld & Melnick|[T99TI : lHriinphrevsll2007l) . The sub- 
millimetre water maser transitions within the ground-vibrational 
state have been detected towards CSEs of O-rich late-type stars 
as well as from low- and high-mass SFRs, either tracing shocked 
regions or arising from the steadilv expanding CSEs (e.g. Ivison 
et al. 1998; M elnick et al.lll993h . The different water transitions 
observable in the ALMA frequency range within the ground- 
vibration al state and the V2 = 1 bending mode are listed in the 
Table |2] (Picke tt et alJ fl998). The rotational maser transitions 
within the vibrationally excited level V2 = 1 are thought to arise 
from regions where the physical conditions are similar to those 
that invert the level popu lation of the SiO rotational transitions 
dAlcolea&Mentenlll993h . 

In contrast to the case of SiO, there are not many references in 
the literature reporting the detection of H2O maser polarized ra- 
diation in the ALMA frequency range. To date, onlv lHarwit et alj 
(2010) have succeeded in measuring linear polarization of water 
maser radiation at 620.7 GHz. Unfortunately, the very low atmo- 
spheric transmission around 621 GHz prevents the detection of 
this maser transition from ground-based telescopes. 
The excitation conditions of the submillimetre water masers 
seem to be a subset of the broad excitation conditions gen- 
erating the 22.2 GHz maser line dNeufeld & Melnickl [1991: 
[Humphreys 2007). The dominating pumping process depends 
on the characteristics of the region where the maser emission 
arises; either from post-shock regions, where the inversion of 
the population of the molecular states is mainly collisional, or 
from the steadily expanding CSEs, where the pumping process 
is mainly thought to be caused by infrared photons of warm 
dust emission. Observations of 22.2 GHz H2O maser polarized 
emission have been used to probe the magnetic field strengths 
withi n the steadily expandin g CSEs of AGB and supergiant stars 
(e.g. I Vlemmings et al.ll2005l) . In addition, polarized water maser 
emission has been detected towards the high-velocity outflows of 
the so-called water-fountains, with spectral features displaying 
high levels of fractional linear polarization (pi > 5%) and un- 
usually broad velocity ranges. These observations have probed 
magnetically collimated "jets" that appear during the post- AGB 
phase, and are thought to be the precursors of bi polar (multipo- 
lar) planetary nebulae (Vlemmings et al. 2006a; Perez-Sanchez 
etal. l2011h . 

3.3. HCN 

Strong maser emission has also been detected from HCN in the 
CSEs of several carb on-rich (C-rich ) AGB stars in the ALMA 
frequency range (e.g. lBieging||200Tl; ISchilke & Mentenll2003l) . 
Table [2] summarizes the different HCN transitions that can be 
observed as maser emission in the submillimetre wavelength 
regime dMiiller et alj 1200 lb . Linearly polarized HCN maser 
emission of the 89.087 GHz, 7=1-0 transition, within the 
(0,2°,0) vibrationally excited state, has been detected at approx- 
imately 20% towar ds the innermost re gion of the CSE of the 
C-rich star CIT 6 (iGoldsmith et al.l fl988>. The masing region is 
thought to be located between the photosphere and the inner ra- 
dius of the expanding envelope, similar to the SiO masers in the 
oxygen-rich late-type stars. The pumping of the HCN 7=1-0 



Table 3. Fixed parameter values assumed in our numerical mod- 
els. 
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Notes. <n) Values calculated using Equation [2] for a magnetic field 
B = 1 G. {b) Magnetic field of the same order of magnitude as the values 
detected towards SFRs and CSE of late-type stars for regions at densi- 
ties of n Hl = 10'° ctrr 3 in the case of SiO and HCN, and n Hl = 10 s cm" 3 
in the case of H2O. 



transition is more likely caused by the absorption of infrared 
photons, not by collisional processes, though a combination o f 
both processes cannot be ruled out (Goldsmit h et al.l fl988). 
Anisotropics in the population inversion of the masing transi- 
tions might affect the fractional linear polarization, but more ob- 
servations are needed to determine the exact cause of the polar- 
ization. 



4. Model results 

Maser emission can be affected by several non-Zeeman pro- 
cesses that can enhance or even produce linear and/or circular 
polarization, such as anisotropic pumping or the change of the 
quant ization axis along the amplification path ( Wiebe & Watson 
1998). Therefore, it is necessary to solve the radiative transfer 
equation for a particular rotational transition to leam whether the 
polarization detected can be correlated to a large-scale magnetic 
field permeating the masing region. Furthermore, it is necessary 
to determine if the involved rotational transition can produce ap- 
preciable levels of pl when the quantization axis is defined by 
the magnetic field direction, i.e. when the Zeeman frequency gQ 
is higher than the stimulated emission rate R. The inequalities 
gQ > r and gQ, > R allow us to analyse the polarization ob- 
served injenns ofthe ordinary population of the magnetic sub- 
states (IWatsonl l2008). In the case of SiO the decay rate F listed 
in Table [3] corresponds to the rate for the radiative decay from 
the first vibrationally excited s tate to the ground-vibrational state 
(Nedoluha & Watson 1990a). For water, this value is roughly 
the inverse lifetime for infrared transitions of th e 22 GHz, H2O 
maser transition dNedoluha & Watsonl [l990b). 

In the case of 

HCN, this value corresponds to the decay rate associated to a 
pumping process dominated bv infrared radiation (Goldreich & 
Kwan ll974l) . Considering the T and sQ values in Table [3] it is 
clear that for the case of interest the criterium gQ > F is satis- 
fied. To evaluate whether the stimulated emission rate satisfies 
the condition gQ > R for a particular transition when pi is suf- 
ficiently large to b e detected, we have used a radia tive transfer 
code a dapted from lNedoluha & Watsonl dl992h and IVlerrim ings 
J2002T) . The stimulated emission rate R is given by 



R 



AkT h AQ 
Anhv ' 



(1) 



where A is the Einstein coefficient of the involved transi- 
tion, which is calculated and listed for most of the transitions in 
Tables[T]and[2j k and h are the Boltzmann and Planck constant, v 
is the maser frequency, 7), is the brightness temperature and AQ 
is the relation between the real angular size of the of the mas- 
ing cloud and the observed angular size. Although the quantities 
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Tb and AQ are related to the observed intensity of the maser, 
it is a difficult task to constrain their value directly from obser- 
vations. Brightness temperatures of up to ~ 10 15 K have been 
measured for water masers in SFRs, whereas estimated values 
for AO are ~ 10~' - 10~ 2 s r for water masers de tected towards 
the CSEs of late-type stars (Richards e t al.ll201 It Vlemmings & 
van Langevelde l2005h . and in some cases AQ ~ 1(T 5 sr for the 
masers detected towards SFRs (Nedoluha & Watsonl [l991l) . 
The fractional linear polarization increases with r^AQ in the 
presence of a magnetic field as long as 9 + 6 cr . Since R increases 
faster than pi when T^AQ increases, the polarization vector of 
linearly polarized maser radiation detected with high brightness 
temperature does not necessarily satisfy the criterium gQ > R, 
and thus cannot always be directly correlated with the direction 
of the magnetic field lines, unless the emitted radiation has a high 
degree of beaming. On the other hand, low brightness tempera- 
tures might result in very low or undetectable values of fractional 
linear polarization. Furthermore, the value of 7/,AQ scales with 
the ratio R/F, and higher values of the molecular decay rate F 
imply a lower level of fractional linear polarization for the same 
T^AQ value. This reflects the dependence of the fractional linear 
polarization level of a particular rotational transition on the sat- 
uration level of the maser radiation. 

The Zeeman frequency gQ determines the energy-splitting of the 
magnetic sub-levels. In a masing region permeated by a con- 
stant magnetic field parallel to the z-axis of the coordinate sys- 
tem, the energy -shifting of the magnetic sub-levels is given by 
hgQm/2 - gp^Bm, where p^ and m are the nuclear magneton 
and the quantum number of the magnetic substate (Nedoluha & 
Watson ll990air Hence the values for gQ. can be estimated using 
the relation 



gQ. 



2gn N B[G] 

n lxio 4 ' 



(2) 



where the factor 1 x 10 4 originates from the conversion be- 
tween the units of Tesla and Gauss. In general, the molecular 
Lande factor (g) is different for each magnetic sub-level. The 
data in the literature are very limited, and there are no reported 
g-factor values for all the different rotational transitions of HCN 
nor H2O. The Lande g-factor for SiO has minor differences (less 
than 1 %) for the v = 0, v = 1 and v = 2 vibrational transi- 
tions dLandolt-Bomsteinll 19821) . The assumed ^-factors for the 
diffe rent molecules are list ed in Table [3] In the case of H2O, 
iNedoluha & Watsonl d!992h calculated the gQ values consider- 
ing hyperfine splitting for the 616 and 523 rotational states of wa- 
ter. Nevertheless, in the present case, we have assume the lower 
limit suggested for g as it would be without hyperfine splitting 
(iKukolichl 1 969h . to calculate gQ for the different H2O rotational 
transitions that we modelled. For HCN we assumed the lower 
g factor value reported for the (0,1 lf ,0) vibrational state, which 
correspo nds to a magnetic field parallel to the molecular symme- 
try axis (Goldsmith et al. 1988 and references therein). Because 
for closed-shell, non-paramagnetic molecules the response to a 
magnetic field is weak, we consider the g-factor values listed in 
Table [3] as a conservative choice to constrain a minimum value 
for gQ for the different molecular species. The gQ values pre- 
sented in Table [3] correspond to B= 1 G in Equation [2] and can 
be scaled according to magnetic field strengths reported for the 
different masing regions in both SFRs and CSEs of (post-)AGB 
stars. 

Single-dish observations have revealed average magnetic field 
strengths of ~ 3.5 G in the SiO maser region, whereas for H2O 
the measured values range between 100-300 mG for the CSEs 



Table 4. Model results: Values of pi assuming 9 = 90° andgQ = 
107? for rotational transitions of the v = 1 vibrational state of SiO 
and the ground-vibrational level of H2O. For HCN, the 7=1-0 
belongs to the vibrationally excited V2 = 2° state, whereas the 
higher 7 transitions listed are from the V2 = 1 c vibrationally 
excited level 



SiO 






H,0 




HCN 




Ju - Jd 


Pl 


Ju 


-h 


Pl 


Ju - Jd 


Pl 


1-0 


-0.32 


3,3 


- 2-20 


-0.19 


1-0 


-0.33 


2-1 


-0.23 


5l5 


-4 22 


-0.16 


2-1 


-0.27 


3-2 


-0.18 


643 


- 550 


-0.15 


3-2 


-0.22 


4-3 


-0.15 








4-3 


-0.19 


5-4 


-0.14 












6-5 


-0.13 













of AGB and supergiant stars, and between 15-150 mG at densi- 
ties of = 10 8 - 10 11 cirT 3 in SFRs. For our models we as- 
sumed the values of the magnetic field strength listed in Table[3] 
This allows us to analyse our results as a function of the mini- 
mum gQ values for the three molecular species modelled. Thus, 
using Equation|2]for the corresponding magnetic field in Table[3] 
we calculate the gQ in order to determine whether the inequality 
gQ > R is satisfied or not. 

5. Analysis 

5.1. SiO maser 

We present the results of our models for the SiO rotational tran- 
sitions 7 = 5-4 and J — 2 - 1 of the v = 1 vibrational state 
in Figure [T] The results are for four different 9 values and the 
vertical line corresponds to a T/,AQ where gQ = IQR. The ideal 
case of 9 — 90° determines the maximum fractional linear po- 
larization value when a magnetic field permeates the masing re- 
gion. The maximum fractional linear polarization for rotational 
transitions of the first SiO vibrationally excited state are listed 
in Table [4] These upper limits were established using the corre- 
sponding gQ value listed in Table [3]for each molecular species, 
and subsequently finding the T^AQ values where gQ = IQR. 
According to our results, it is possible to generate pi values of 
up to 13% for the 7 = 5-4 rotational transition, and of up to 
23% for the 7 = 2-1 rotational transition without consideration 
of non-Zeeman effects (Figure [TJ, while the much higher values 
of pi detected for the SiO masers require anisotropic pumping. 
Therefore, although anisotropic pumping has a strong impact on 
the fractional linear polarization level, submillimetre SiO maser 
features observed with brightness temperatures < 10 9 - 10 10 K 
fulfil gQ > R and could in principle be used to trace the structure 
of the magnetic field permeating the SiO masing region. 

5.2. H 2 maser 

To determine the level of fractional linear polarization that water 
maser transitions can reach when the masing region is perme- 
ated by a large-scale magnetic field, we ran numerical models 
of the rotational transitions within the ground-vibrational state, 
using the assumptions described in section IBTl together with the 
corresponding parameters in table [3] The results are presented 
in Table |4] In Fig. [2] we present the results for the 183 GHz and 
325 GHz lines for four different 9 values. Our results suggest up- 
per limits for fractional linear polarization of up to 19% and 16% 
for the 183 GHz and 325 GHz lines, respectively. Therefore, 
considering the values we assumed for the input parameters 
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215.59 GHz J = 5-4 SiO v=l 183 GHz 3 13 -2 20 H 2 
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T b Afi (K sr) T b AQ (K sr) 




Fig. 1. Fractional linear polarization as a function of TbAQ for 
the 215 GHz, 7 = 5-4 (top), and the 86.2 GHz, 7 = 2-1 (bot- 
tom) SiO rotational transitions within the v = 1 vibrational level 
for four values of the angle between the magnetic field lines and 
the direction of propagation of the maser radition, 0, using the 
parameters listed in table[3]for SiO. Positive and negative pi val- 
ues mean that the vector of polarization is either perpendicular 
or parallel the magnetic field lines projected on the plane of the 
sky. The vertical line sets the TbAQ value such that gQ = 1QR. 

listed in Table[3] maser features with 10 K< 7/, < 10 11 K can 
produce observable fractional linear polarization levels that can 
be used to determine the magnetic field morphology. However, 
if Tt, is larger, the observed polarization direction will be rotated 
away from the projected magnetic field direction on the plane of 
the sky. 



Fig. 2. Fractional linear polarization as a function of TbAQ. for 
the 183 GHz, 3 i3 -2 20 (top), and the 325 GHz, 5i 5 -4 22 (bottom) 
H 2 rotational transitions, both within the ground- vibrational 
level, for four values of the angle between the magnetic field 
lines and the direction of propagation of the maser radiation, 0. 
We considered the magnetic field strength and decay rate listed 
in table [3] for H 2 0. Positive and negative pi values mean that 
the vector of polarization is either perpendicular or parallel the 
magnetic field lines projected on the plane of the sky. The verti- 
cal line sets the TbAQ value such that gQ = 10/?. 



5.3. HCN maser 

In Table|4]we present our results for the maximum fractional lin- 
ear polarization that the v 2 = 2°, 7=1-0 rotational transition 
and higher 7 transitions within the v 2 = 1 1[ vibrational state can 
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Fig. 3. Fractional linear polarization as a function of 7/,Afi for 
7 — 1 - 0, V] = 2° rotational transition (top), and three HCN 
rotational transitions of the vibrationally excited state vi = 1 
(bottom), assuming a magnetic field oriented perpendicular to 
the direction of propagation of the maser radiation. We consid- 
ered the magnetic field strength and decay rate listed in table 
[3] for HCN. Under these assumptions, the vertical lines set the 
maximum linear polarization that the different maser lines can 
reach satisfying the condition gQ = 10R. 



produce. In Figure [3] we show our results assuming a magnetic 
field oriented perpedicular to the direction of propagation of the 
maser radiation, and the parameters listed in table [3] for HCN. 
Fractional linear polarization of up to ~ 33% can be generated 
in the case of the 7=1-0 rotational transition for gQ = 10R. 
This value decreases for rotational transitions with higher 7 val- 
ues, as expected. The values of Tjfl scale with R/F, and conse- 
quently, lower F values will increase the fractional polarization 
level that a particular rotational transition can reach while still 
satisfying gQ > R. Nevertheless, according to our results, HCN 
maser emission towards C-rich AGB stars can be generated with 
high fractional polarization values within the gQ. > R regime. 
Thus, the vector of polarization detected from HCN maser ob- 
servations could be used to probe the magnetic field structure 
around C-rich AGB stars, even if the level of polarization has 
been enhanced by non-Zeeman effects. 



6. Observing maser polarization with ALMA 

The temperature and density conditions that favour the popula- 
tion inversion of different masing molecules in SFRs and CSEs 
of (post-) AGB stars are not the same. Most problably these con- 
ditions also differ for different rotational transitions. Therefore, 
detecting polarized maser emission from multiple rotational 
transitions and from multiple molecular species with ALMA will 
enable us to trace the magnetic field structure throughout ex- 
tended regions around those sources. The accuracy of the ALMA 
polarimetry will allow us to detect polarized emission of 0. 1 % 
of the detected Stokes I. Therefore, fractional linear polarization 



Table 5. Sensitivity values at the different maser frequencies of 
SiO, H 2 and HCN within the ALMA frequency range. The 
<x listed corresponds to the sensitivity achieved after lh of on- 
source time, with 0. 1 km/s for a source that reaches a maximum 
elevation of 59°. 



SiO 




H 2 




HCN 




V 


cr 


V 


O" 


V 


cr 


(GHz) 


(mJy) 


(GHz) 


(mJy) 


(GHz) 


(mJy) 


86.243 


4.58 


183.310 


38.37 


89.0877 


4.49 


129.363 


4.32 


325.153 


39.58 


177.238 


4.97 


172.481 


4.74 


439.151 


59.25 


267.199 


4.12 


215.596 


3.63 






354.461 


4.99 


258.707 


4.76 











can be detected in short observations with very good signal-to- 
noise ratios. To probe whether the measured polarization of a 
maser traces the magnetic field permeating the masing region, it 
is necessary to calculate the brightness temperature of the spec- 
tral feature and evaluate if the product T b AQ satisfies the in- 
equality gQ. > R for the observed emission. For detected spectral 
features it is possible to constrain T b by considering the equation 



[K] [Jy]l[mas 2 ]/ 



(3) 



where S (v) is the detected flux density, S is the maser angular 
size and f v is a constant factor that includes a prop ortionality 
factor obtained for a Gaussian shape (Bu rns et al 1l979l It scales 
with frequency according to the relation 



£, = 6.1305x 10' 



i— V 

\GHz/ 



-K. 



(4) 



We calculated the rms value at the different maser frequen- 
cies that we modelled using the ALMA sensitivity calculator 
for an array of 50 12-m antennas. The results are presented in 
Tableland correspond to observations of 1 h of on-source time 
with 0.1 km/s of spectral resolution. To estimate the brightness 
temperature from observations that do not resolve individual 
features, it is necessary to assume a value for the size of the 
masing region generating the spectral feature. VLBI studies 
have revealed masers to be very compact spots with typical sizes 
of 1 AU. Based on observed peak-flux density values for maser 
emission of SiO, H2O and HCN, assuming a size of 1 mas, 
i.e. a source distance of 1 kpc, we here give a few examples of 
potential ALMA observations of maser polarization. 

SiO masers: SiO maser emission has been reported dis- 
playing flux densitv values between 7.4 Jv and 64 Jv (Shinnaga 
et al. 120041) for the v = 1,7 = 5-4 rotational transition. 
Moreover, Barvainis & Predmorel (119851) reported flux density 
values between 210 Jy and 625 Jy for the SiO 7 = 2-1 
rotational transition of the same vibrationally excited state. As 
an example, we investigate the detection of a maser spectral 
feature with a peak- flux density of 10 Jy of the SiO 7 = 5-4 
(v = 1) rotational transition. For these values, Equation [3] 
gives T b = 1.32 x 10 8 K. Comparing this with the T b AQ 
value, which corresponds to the limit of gQ - 10R (Figure 
[TJ, maser emission with beaming angles AQ < 2.84 x 10 _1 sr 
will have a linear polarization that can be used to trace the 
magnetic field direction. But usually SiO maser spectra exhibit 
blended components. If this is the case, the flux density ob- 
served corresponds to a number of maser spots with similar 
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line-of-sight velocities increasing the observed S(v). Therefore, 
the brightness temperature derived using Equation [3] could be 
overstimated, but still could be used to set an upper limit value 
for T/,. To reach the regime where gQ > R, it is necessary to 
associate a AO value to the brightness temperature derived from 
the S(v) of the blended components. Hence, since individual 
maser spectral features might have lower T/, than the blended 
feature which contains it, even unbeamed maser radiation could 
still place the linear polarization detected within the regime 
where gQ > R is satisfied. 

Finally, according to the ALMA sensitivity calculator, a 5<x 
detection of 1% of fractional linear polarization of SiO maser 
emission with low peak-flux densities (few tens of Jy/beam) 
requires short on-source observation time (t < 1 h). 

H2O masers: Submillimetre H2O lines have been detected 
with flux densities from a few tens to several thousands of 
Jy. For the minimun peak- flux density (243 Jy /beam) reported 
for the 183 GHz by Ivan Kempen et all (l2009h . the brightness 
temperature we obtain using Equation [3] is Tt, = 4.43 x 10 9 K. 
Hence, H2O maser emission detected with beaming solid 
angles AQ > 1.6 x 10~ 3 sr can generate fractional linear 
polarization levels higher than 1% (Figure |2), satisfying the 
condition gQ > R. The 1% of linear polarization level of 
the weakest 183 GHz maser feature detected bv van Kempen 
et al. d2009l) can be easily detected in very short integration 
times. Unfortunately, many H2O maser lines are affected by the 
atmospheric PWV. For maser spectral features with peak-flux 
densities of the order of 20 Jy/beam at 439.151 GHz, for 
instance, a 3cr detection of 1 % will require ~ 1 h of on-source 
observation time. 

HCN masers: HCN / = 1 - 0, (V2 = 2°) maser emission 
at 89.0877 GHz has been detected with fractional linear polar- 
ization of 20% for a spec tral feature wit h flux density of 38 Jy 
(Goldsmith_et_aL| [l988). furthermore, Lucas & C ernicharo 
(1989) leported the detection of HCN V2 = 1 c maser emission 
with a peak-flux density ~ 400 Jy at 177 GHz. ISchilke et ail 
(2000) reported the detection of the vibrationally excited HCN 
(04°0) 7 = 9-8 maser line near 805 GHz with a flux density of 
-1500 Jy. 

For emission reported by Golds mith et al.l (119881) . assuming a 
magnetic field perpendicular to the propagation of the maser 
radiation, 20% of fractional linear polarization corresponds 
to T h AQ = 7.65 x 10 5 K sr (Figure QJ top), which places the 
polarization detected in the regime where gQ > R. On the 
other hand, considering a spectral feature of S(v) = 38 Jy, 
Equation [3] gives 7j = 2.94 x 10 9 K, implying a beaming angle 
of 2.6 x 10~ 4 sr, assuming a maser of 1 mas. A beaming angle 
AO. ~ 10 -2 sr is needed for this maser to satisfy gQ. > IQR and 
be a realiable tracer of the magnetic field within the masing 
region. However, if the emission consists of a contribution from 
multiple maser spots, the brightness temperature of the spectral 
feature will be overstimated, as discussed previously for blended 
SiO maser lines. Finally, as in the case of SiO maser emission, 
1 % of linear polarization could be observed with very low rms 
values in short on-source observation times (Table |5). 

7. Conclusions 

We have run numerical models to calculate the fractional linear 
polarization pi of maser emission generated by the interaction 
of a magnetic field with the different rotational transitions of 
SiO, H 2 and HCN within the ALMA frequency range. The 



fractional linear polarization was calculated as a function of 
ThAQ., a quantity that can be related to the stimulated emission 
rate R of the involved transition. Considering both the minimun 
value of the Lande g-factor for each molecular species and 
a suitable magnetic field strength for the different masing 
regions, we found the maximum pi that the analysed rotational 
transitions can reach while satisfying both conditions gQ > R 
and gQ > T. Meeting these criteria allow us to use the detected 
vector of polarization as a tracer of the magnetic field in the 
masing region, even if the polarization observed has been 
affected by non-Zeeman effects. 

According to our results, SiO, H2O and HCN submillimetre 
maser emission can be detected with observable fractional 
polarization levels (> 1%) in the regime gQ > R. But es- 
pecially for the strongest masers, a careful analysis of the 
brightness temperature is needed to confirm that the maser 
polarization is still in this regime. Thus, observing with ALMA 
full-polarization capabilities will enable us to use polarized 
maser emission as tracer of the magnetic field structure towards 
SFRs and CSEs of (post-)AGB stars. Depending on the maser 
spectral features detected, both the brightness temperature and 
the beaming solid angle can be better constrained, leading to a 
more accurate determination of the direction of the magnetic 
field in the masing region. 
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